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We perform a wavelet analysis of the temperature and polarization maps of the Cosmic Microwave 
Background (CMB) delivered by the WMAP experiment in search for a parity violating signal. Such 
a signal could be seeded by new physics beyond the standard model, for which the Lorentz and CPT 
symmetries may not hold. Under these circumstances, the linear polarization direction of a CMB 
photon may get rotated during its cosmological journey, a phenomenon also called cosmological 
birefringence. Recently, Feng et al. have analyzed a subset the WMAP and BOOMERanG 2003 
angular power spectra of the CMB, deriving a constraint that mildly favors a non zero rotation. By 
using wavelet transforms we set a tighter limit on the CMB photon rotation angle Aa = —2.5 ± 3.0 
(Aa = —2.5 ± 6.0) at the one (two) o level, consistent with a null detection. 



I. INTRODUCTION 



The CMB is one of the primary experimental windows to the early universe. Recent observations have reached 
remarkable precision. When combined with other complementary cosmological datasets, the WMAP three year 
(hereafter, WMAP3) observations [l[ convincingly support the so-called standard model of structure formation (2j. 
However the CMB may provide further information. In principle, one may use the background photons to constrain 
new physics beyond "standard" models. A positive answer might be provided by the study of CMB polarization 
(CMBP), whose observations currently mark the experimental frontier of the held. Pioneering observations, including 
DASI [|, CBI Q, BOOMERanG 2003 [| (hereafter B03), MAXIPOL Q, and WMAP itself Q have yielded detections 
of the CMBP over a wide range of angular scales. Within the next decade, ground or space-based experiments may 
detect via the CMBP a signal from primordial gravitational waves Q , thus constraining the energy scale of the inflation 
and probing particle physics well beyond the capability of any conceivable terrestrial accelerator. The CMBP can 
also provide information on symmetry-violating physics beyond the Lee- Yang parity (V) breaking that is central to 
the standard model, yet not observable through CMB anisotropics due to their charge blind character. In general, 
the breakdown of spacetime symmetries is a potential tracer of new physics [9(. Several models exist that predict 
non-standard V and CP violations ('C standing for charge conjugation), as well as CPT violations ('T' being time 
reversal) and the related (through the anti-CPT theorem (To| ) breakdown of Lorentz invariance. A number of tests 
have been suggested and (in many cases) performed, either in terrestrial and orbital laboratories (ll| or through 
cosmological observations [HI, [l3[ . These violations might also have a measurable imprint on the observed CMBP 
pattern, whose statistical properties are constrained by the assumption of symmetry conservation. 

For a sky direction n, a polarized map of the CMB is usually given in terms of total intensity (or temperature) T(n) 
and linear polarization Stokes parameters Q(n) and U(n). The T field can be decomposed into scalar (S) spherical 
harmonics Yi m (n), obtaining the coefficients af m - Q and U are components of a symmetric, trace- free rank 2 tensor, 
and are expanded in tensor spherical harmonics Y^ n {n) and Y^ l {n) with coefficients af m and af m , respectively. These 
correspond to scalar (gradient-like) 'G' and pseudo-scalar (curl-like) 'C modes [8]. Under hypothesis of Gaussianity 
and isotropy, the statistical properties the CMB are described by two point correlation functions on the sphere, 
whose Legendre transforms define six angular power spectra: Cf z = {af m {af m )*) with Z, Z' = {T,G,C}. If the 
physics controlling CMB fluctuations is parity conserving the cross spectra Cj c and Cf c must vanish due to the 
different handedness of the C and (S,G) harmonics. Therefore, if the standard cosmological model holds ,we should 
expect no relevant information from TC and GC . On the other hand, detection of non-zero primordial TC and/or 
GC may probe fundamental physics in the early universe, such as the presence of a primordial homogeneous [l4| or 
helical [151 ] magnetic field which would induce Faraday rotation and non-zero TC correlations. Parity-asymmetric 
gravity dynamics during inflation may generate a discrepancy among left and right-handed gravitational waves, so 
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that TC and GC are non-zero [l6j]. Particle physics models with non-standard parity- violating interactions also 
predict non- vanishing TC and CG signals (l7| . 

In this paper we focus on a class of models that exhibit parity violations in the photon sector. A Chern-Simons 
term is introduced in the effective Lagrangian (l3j : 

where F^ v is the Maxwell tensor and the 4-potential. The 4-vector p M may be interpreted as the derivative of 
the quintessence field or the gradient of a function of the Ricci scalar [l8[ • In either case a V violation always arises 
provided that po is non-zero, while C and T remain intact. Hence, CV and CPT symmetries are also violated, as well 
as Lorentz invariance, since p^ picks up a preferred direction in space-time. The net effect on a propagating photon 
is to rotate its polarization direction by an angle Act, hence the name "cosmological birefringence". Historically, 
the effect has being constrained by measuring polarized light from high redshift radio galaxies and quasars (l3l . [l9| . 
Obviously, the CMB photons would also be affected and, due to their longer journey, may get a larger rotation. A 
consequence for the CMB pattern is the mixing of G and C modes: the TG and GC correlations still vanish at last 
scattering surface, but the observable CMB spectra are distorted as [HI [20j: 
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where the primed quantities are rotated. In [2l[, the TT and TG power spectra measured by WMAP3 together with 
all six spectra measured by B03 have been used to perform a global fit, yielding a mild detection for a non zero 
rotation (but see also [24| for a similar analysis restricted to the B03 power spectra and [25[ where constraints on the 
coupling between the quintessence and the psudoscalar of electromagnetism are derived, based again on B03 data). 
Using the same data set, a similar result has been found in [22| . and used to constrain a specific baryo/leptogenesis 
model, while an interaction between the neutrino asymmetry and a term Chern-Simons term has been proposed in 
[23| as a possible explanation for the result found [2l| . Forecasted constraitns on Aa for high sensitivity experiments 
such as Planck or CMBpol can be found in [26| . 

Here we constrain Aa with a wavelet analysis. A rotation of the photon polarization direction leaves an imprint on 
each resolution element (or pixel) of the Q and U maps, and a map-based estimator appears appropriate. Wavelets 
are a natural choice because they allow for multi-scale pixel analysis. We compute the wavelet cross-correlation 
coefficients for TC and GC to build a goodness of fit estimator that we apply to the WMAP3 {T, Q, U} maps. Our 
analysis is complementary to that of [21(, where the information for TC and GC comes from B03. The two analyses 
differ in the method and (substantially) in the data set (the only overlap being the WMAP3 temperature map). In 
the following, we derive more stringent limits on Aa by adapting the wavelet formalism to tackle polarization, a point 
that has not been addressed to date in the CMB literature. 

The plan of the paper is as follows: in section [TT] we describe our wavelet based method to constrain Aa: section Hill 
we define a suitable estimator and apply it to WMAP3 data, making use of numerical simulations. Finally in section 
IIVI we draw our conclusions. 



II. A WAVELET STATISTIC FOR TEMPERATURE AND POLARIZATION 

Given a position X, wavelets are filter functions ^(X;b, R) that also depend on a characteristic scale R and 
translation b. They provide scale-varying transforms that remain localized in pixel space. Moreover, they consist 
of an infinite set of basis functions, thus providing some freedom of choice in matching their functional form to the 
target signal. Several authors have exploited this flexibility as a powerful tool in CMB data analysis. Wavelets have 



been used for denoising |27|, point source extraction [28(, foreground removal [29( and for detecting the integrated 
Sachs Wolfe effect [3(| ■ Since the wavelet transform preserves linearity, its coefficients can be used to constrain the 
statistics of the field at different scales. In particular, the spherical Mexican hat wavelet (SMHW) has been used to 
flag statistical anomalies in the WMAP data [3l| and to constrain primordial non-Gaussianity [13] (other types of 
wavelets have been shown to be sensitive to yet different anomalies, e.g.(33j). SMHW are generated from ordinary MH 
wavelets through a stereographic projection on the tangent plane [34j that is known to preserve their basic properties 
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The SMHW is defined as: 
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where y = 2 tan 8/2 (8 is the polar angle), R is the scale of convolution, and N(R) a normalization factor. For a T 
map, the wavelets coefficients are: 

W T (R,n) = J dfl'T(n + h')^(e',R) . 
This convolution can be performed in harmonic space: 
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where ^>i(R) are the Legendre expansion coefficients of the SMHW. Handling polarization requires more care, since 
Q and U are not rotationally invariant, being components of the rank 2 tensor P a j, [8j. By taking the covariant 
derivatives of P a b, one can build two quantities that are rotational invariant and hence decomposed by S harmonics. 
This leads, again, to the G and C coefficients: 
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Here ':' stands for the covariant derivatives on the sphere, e is the Levi-Civita trace- free antisymmetric tensor and 
Ni a normalization factor Q. We similarly define SMHW coefficients as: 



W G {n,R) = J dn'P^ b (n + h')^(0',R) 



W c (n, R) 
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Note that we never explicitly compute derivatives, since the integrals can be performed in harmonic space (c.f. [6]), 
provided we divide out the factor Ni. Finally, we consider the pixel-pixel cross correlation of the SMHW coefficients 
as our main statistic: 
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where Z = {T, G} and V is a volume normalization that can be taken to be proportional to the total number of pixels 
N p . The quantities in cq. [7] possess the same P symmetry of the usual harmonic cross spectra: they can be non-zero 
only if parity conservation is violated. 



III. NUMERICAL SIMULATIONS AND RESULTS 

To constrain Aa, the following scheme was employed. We modified the Healpix package [36[ to generate a set of 
Monte Carlo (MC) simulations for {T,Q,U} maps containing a CMB signal whose polarization pattern is rotated 
according to eqs. [I][5j we use the WMAP3 best fit model as the (unrotated) fiducial angular power spectrum. The 
signal maps were smoothed according to the WMAP3 optical transfer function. We also simulated noise maps 
consistent with the WMAP3 instrumental properties. We add to each signal map a noise realization consistent with 
the WMAP3 instrumental properties. Simulating noise in TQU maps is more complicated than for T only, because 
the noise values of different Stokes parameters within a given pixel are usually correlated. For WMAP3, T is very 
weakly correlated with Q and U, so this coupling can be safely neglected [37j • On the contrary, in order to obtain 
accurate results one has to take into account the correlations between Q and U. The WMAP team has released 2x2 
effective hits arrays (hereafter, iV bs where the off diagonal elements represent the (QU) inter pixel correlation (off 
pixel correlations are very weak and can be neglected); these matrices are given for each differential assembly (DA) 
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FIG. 1: MC means for X T (see text), for —8° < Aa < 8°, step of 1° (dashed lines refer to positive Aa, the middle line 
is for Aa = 0, i.e. V is conserved). The scale R is given in arcminutes. The shaded region shows the la range for Aa = 0. 
Experimental points (WMAP3) are shown as crosses. 



and for each observation year. To simulate noise maps, the following scheme is employed: for a given DA and for each 
pixel, we add the N Q b s arrays for different years. The three year, noise maps for each DA i are then simulated as: 



Qi. 

u. 
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where p identifies a given pixel, (tqu is the nominal DA polarization sensitivity, as provided by the WMAP team [37| 

1/2 

and N obs is the Choleski factor of N bs- The resulting noise plus signal maps for each DA are then weighted averaged 
to form a combined map comprising all DA's in the Q,V and W bands: 
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The combined map for T is computed by using the standard, scalar version of the procedure above: 
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To minimize residual foreground contamination we chose to use a rather conservative mask, the intersection of the 
KpO and P02 sky cuts [J. The masked maps are downgraded in resolution to 13.6'. We then compute the wavelet 
coefficients W z (p,R) (where Z = {T,G,C}) over the discretized sphere. We consider 17 wavelet scales from 14' to 
100'. 1 To avoid boundary effects, we widen the sky map up to a fraction 2.5 R [3lT l32l ]. Finally, we consider the 
goodness of fit statistics x 2 (Aa) = y T C^ a Y, where Y = X WMAP -X(Aa). X WMAP is computed over the Q+V+W 
foreground-cleaned, optimally-weighted data map, while the mean (barred) quantities are derived from ~ 2000 MC 
simulations. The covariance matrix C^a is estimated over a fresh set of simulations (~ 4000). 

In fig. Q] we show, as a function of R, the MC mean values for X TC , for —8° < Aa < 8° with a step of 1°. The 
crosses are experimental points from WMAP3 and the shaded region is the la range, centered about the Aa = case. 
In fig. [2] we show the same for X GC , with the means in the range —16° < Aa < 16° with a step of 2°. As expected for 
WMAP3, TC has a significantly larger signal to noise than GC. The means are computed over noisy simulations, but 
closely reproduce the ensemble predictions that can be derived from eq. [7J To show that our estimator is unbiased, we 
simulated a further MC set with given Aa and checked that the means of the x 2 estimates reproduce the input values 
with high accuracy. Throughout our analysis, we keep the dependence of Aa in the estimator's covariance matrix 
(but find no significant change in our results if we drop this dependence: for WMAP3, our estimator's covariance is 
completely dominated by noise). 

In fig. [3] we show the likelihoods of WMAP3 data for TC and GC. GC contributes very weakly to the joint 
likelihood C oc exp(— % 2 /2). We estimate Aa = —2.5 ± 3.0 and Aa = —2.5 ± 6.0 at la and 2a confidence limit 
respectively. Thus, we find no evidence of parity violation from the WMAP3 data. These limits are slightly tighter 
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FIG. 2: Same as fig. [T] but for X GC and -16° < Aq< 16°, step of 2°. Note the lower signal to noise. 
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FIG. 3: Likelihood functions of the cosmological birefringence angle Aa for CG (dotted) and TC (dashed), computed from the 
wavelet estimator on WMAP3 data. The solid (blue) line is the global, covariance weighted, likelihood 

than those given in (2l| , where a marginal detection for a non-zero Aa is claimed and seems to be driven by the GC 
B03 data. To show that our conclusions do not depend on the particular fiducial power spectrum chosen (provided 
it is reasonable), we have allowed the latter to vary between the ±lcr experimental limits set by WMAP3, finding 
fully consistent results (no detection, very similar limits). This procedure extends to polarization the test suggested 
by |31| for temperature data. 

As a futher consistency check, we compared our estimator with a similar one, built using the angular power spectrum 
rather than wavelets. In the case of pure signal, under the assumption of Gaussianity and statistical isotropy of the 
observed field, we expect the two approaches to provide similar constrains. To show this is indeed the case, we have 
repeated the procedure of section HTT1 performing a Monte Carlo simulation over 1000 realizations in the ideal case of 
pure signal. The % 2 using the cross spectrum Cj c can be calculated analytically by: 

X \Aa)=Y J (Ct-Crf/^ (11) 

where as usual the prime identified rotated spectra, and the cosmic variance is given by Q : 

^ = [{cr? + crci cc \/2i+i 

and the maximum multipole in the sum is l m ax — 500 roughly consistent with the maximum resolution emploied in 
the wavelet analysis. In figure|4]we show \ 2 a $ a function of Aa against the null hypothesis for wavelets and the cross 
spectrum Of . The two methods give very similar results, as expected in this ideal case. 

IV. CONCLUSIONS 

In summary, we have presented the first application of wavelets to polarized CMB maps, and used it to constrain 
the rotation angle of CMB photons in search of a signature due to cosmological birefringence, an effect connected to 
fundamental symmetry-breaking physics. We find no evidence of such a rotation and present the best upper limits 
to date on CMB data. This result should be compared with [2l|, where a marginal detection for a non-zero Aa is 
claimed. The latter result is mostly based on the B03 data and only makes use of a subset of the WMAP3 dataset, 
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(dotted) and the angular power spectrum Cf c (solid) 

not including the TC correlations from which our results are essentially derived. While WMAP3 has lower signal to 
noise per pixel than B03, the analysis presented here uses data from ~ 60% of the whole sky while the limited useful 
sky coverage of B03 (< 1%) severely limits the statistical power of TC, so the detection in [21| appears to be driven 
from the much harder to measure (and prone to systematic effects) GC correlations 2 . Given the quantity and quality 
of the CMB data anticipated over the next few years, our approach demonstrates that substantially stronger limits 
on parity violation should be feasible. 

V. ACKNOWLEDGMENTS 

We thank V. Antonuccio, O. Dore, M. Lattanzi and A. Palazzo for discussions. This research used resources of the 
NERSC, which is supported by the Office of Science of the U.S. DoE under Contract No. DE-AC02-05CH11231, and 
of the CASPUR supercomputing centre (Rome, Italy). Some of the results presented here have been derived using 
the Healpix package [36[ . 



[1] 
[2] 
[3] 
[4] 

[S. 
[6 

[7] 

[8] 

[9 
[10 

[11 
[12 

[13 



Available from http://lambda.gsfc.nasa.gov 



D.N. Spergel et al, Astrophys. J. in press, [astro-p h/0603449 
J. Kovac et al, Nature, 420, 772, (2002) 
A.C.S. Readhead et al, Science, 306, 836 (2004) 

F. Piacentini et al, Astrophys. J., 647, 833, (2006); T.E. Montroy et al, Astrophys. J. 647, 813, (2006) 



B.R. Johnson et al, |astro-ph/0611394 (2006) 

L. Page, et al, Astrophys. J. in press, astro-ph/0603450 

M. Kamionkowski, A. Kosowsky and A. Stebbins, Phys. Rev. D 55, 7368 (1997); U. Seljak and M. Zaldarriaga, Phys. Rev. 

Lett. 78, 2054, (1997J 

R. Lenhnert, |hep-ph/0611177| 

O.W. Gree nberg, Phys. Rev . Lett. 89, 231602 (2002); Found, of Phys., 36, 1535, (2006) ( |hep-ph/0309309| 
R. Bluhm, |hep-ph/0112318| M. Mewes, |hep-ph/ 0307161 1 
G. Amelino-Camelia et al, Nature, 393, 763 (1998) 

S.M. Carroll, G.B. Field & R. Jackiw, Phys. Rev. D 41, 1231, (1990); S.M. Carroll & G.B. Field, Phys. Rev. D 43, 3789, 
(1991) 

[14] E.S. Scannapieco and P.G. Ferreira, Phys. Rev. D 56, 7493 (1997) 

[15] L. Pogosian, T. Vachaspati and S. Winitzki, Phys. Rev. D 65, 083502 (2002); M. Giovannini, Phys. Rev. D 71, 021301, 
(2005) 

[16] A. Lue, L. Wang & M. Kamionkowski, Phys. Rev. Lett., 83, 1506, (1999); S. Saito, K. Ichiki & A. Taruva, larXiv:0705.370T1 
2007 



2 



Recently, Xia et al. _38| have extended the analysis in 12 II including in their analysis the previously left aside WMAP3 full power 
spectrum dataset. They find Aa = —6.2 ± 3.8deg(l<r), thus confirming a mild detection of a non-zero rotation angle. It would be very 
interesting to understand to what extent this is driven by the B03 dataset. An extension of our analysis to the B03 TQU maps (not yet 
public at the time of writing) is under study. 



7 



[17] D. Maity, P. Majumdar, S. SenGupta, JCAP, 0406, 005, (2004); N.F. Lepora, |gr^qc/9812077[ 1998; K.R.S. Balaji, R.H. 
Brandenberger, D.A. Easson, JCAP, 0312, 008, (2003) 

H. Davoudiasl et al, Phys. Rev. Lett., 93, 201301, (2004); H. Li, M. Li, X. Zhang, Phys. Rev. D 70, 047302, (2004) 
S. M. Carroll and G. B. Field, Phys. Rev. Lett., 79, 2394, (1997); S. M. Carroll, Phys. Rev. Lett., 81, 3067, (1998); B. 
Nodland, J P. Ralston, Phys. Rev. Lett., 78, 3043, (1997); D.J. Eisenstein, E.F. Bunn, Phys. Rev. Lett. 79, 1957, (1997); 
J. P. Leahy, |astro-ph/970428"5| 1997 

B. Feng, H. Li, M. Li, and X. Zhang, Phys.Lett. B 620, 27 (2005) 

B. Feng et a l, Phys. Rev. Let t. 96, 221302, (2006) 
M. Li et a/., |hep-ph/0611192| (2006) 

C. Q. Ceng, S. H. Ho, J. N. Ng, arXiv:0706.0080, (2007) 
V. A. Kostelecky & M. Mewes, |astro-ph/0702379| 

G. C. Liu, S. Lee and K. W. Ng, Phys. Rev. Lett., 97 161 303, (2006) 
J.Q. Xia, H. Li , G.B. Zhao, X. Zhang. larXTv:0708.1111l f2007) 
J. L. Sanz et al, Mon. Not. R. Astr. Soc, 309, 672, (1999) 

L. Cayon et al, Mon. Not. R. Astr. Soc., 315, 757, (2000); L. Tenorio et al, Mon. Not. R. Astr. Soc., 310, 823, (1999); 
J. Gonzalez-Nuevo et al, Mon. Not. R. Astr. Soc. 369 1603, (2006) 
F.K. Hansen et al, Astrophys. J., 648 784, (2006) 

D. Pietrobon, A. Balbi, D. Marinucci Phys. Rev. D 74 043524, (2006) 
P. Vielva et al, Astrophys. J., 609 22, (2004); M. Cruz et al, Astrophys.J., 655, 11, (2007) 

P. Mukherjee & Y. Wang, Astrophys.J. 613, 51, (2004); P. Cabella et al, Mon. Not. R. Astr. Soc, 358, 684, 2005 
J.D. McEwen et al, Mon. Not. R. Astr. Soc. Lett. 371, L50, (2006) 

E. Martinez-Gonzalez et al, Mon. Not. R. Astr. Soc, 336, 22, (2002); L. Cayon et al, Mon. Not. R. Astr. Soc, 326, 1243, 
(2001) 

J. P. Antoine & P. Vandergheinst, J. Math. Phys., 39, 3987, (1998) 

K.M. Gorski et al, Astrophys. J., 622, 759, (2005); web page: http://healpix.jpl.nasa.gov 
N. Jarosik et al, Astrophys. J. in press, astro- ph/0603452 (2006) 
J. Q. Xia , H. li , X. Wang, X. Zhang. larXivT 0710.3325vl. (2007) 



